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Abstract 

Distribution of the two phonon 7 vibrational collectivity in the rotating triaxial odd- A nucleus, 
103 Nb, that is one of the three nuclides for which experimental data were reported recently, is 
calculated in the framework of the particle vibration coupling model based on the cranked shell 
model plus random phase approximation. This framework was previously utilized for analyses of 
the zero and one phonon bands in other mass region and is applied to the two phonon band for 
the first time. In the present calculation, three sequences of two phonon bands share collectivity 
almost equally at finite rotation whereas the K = $7 + 4 state is the purest at zero rotation. 
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I. INTRODUCTION 



One of the properties most specific to the finite quantum many body system, the atomic 
nucleus, is to exhibit both single particle and collective modes of excitation in similar energy 
scale. 

Thanks to recent progress in computer power, exact diagonalization of given effective 
interaction in a huge model space is becoming available for up to medium mass nuclei in- 
cluding those far from stability. But in order to extract physical picture of their dynamics 
from huge numerical information, it is necessary to have recourse to the concept of collec- 
tivity. Among the various kinds of collective modes of excitation, one of those that have 
been well studied and known to exist prevailingly in the nuclear chart is the 7 vibration. 
Vibrational excitations in many fermion systems, in other terms, phonons, are made up by 
coherent superposition of particle-hole or two quasiparticle excitations. Thus, on the one 
hand, multiple excitations are expected if they are really collective. On the other hand, their 
excitation spectra, in particular their anharmonicity reflect underlying nuclear structure. 

In the nuclear chart, collective vibrations in the rare earth nuclei with the mass A =160 
-170 are the best studied. The first experimental information about the two phonon 7 
vibration, denoted by the 27 hereafter, in 168 Er was reported by Davidson et al. [lj based on 
a high resolution 7 ray study following neutron capture. The first theoretical analysis was 
done by Warner et al. I2j in terms of the interacting boson model with s and d bosons but it 

n 

was critically assessed by Bohr and Mottelson |3| . They argued within the general framework 
applicable to deformed nuclei with rotational spectra [l]. Dumitrescu and Hamamoto ({5| elu- 
cidated this problem by means of a macroscopic and microscopic analysis. On the contrary, 
Soloviev and Shirikova [6] argued that collective two phonon excitations do not exist because 
of the strong Pauli principle among nucleons that constitute vibrational excitations based 

n n 

on the quasiparticle phonon model. Matsuo [7| and Matsuo and Matsuyanagi [8] applied the 
selfconsistent collective coordinate method that provides quantized collective Hamiltonian 
starting from the microscopic random phase approximation (RPA) and obtained the result 
that the K = 4 27 state in the triaxially deformed potential exists at the energy about 
2.7 times that of the I7 state. Pipenbring and Jammari [9] also obtained a similar result 
by means of the multiphonon method based on the Tamm-Dancoff approximation. The 
description in terms of the interacting boson model was improved by Yoshinaga et al. loj 
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by introducing the g boson. The definite experimental evidence, the absolute B(E2) value 
between the K = 4 27 candidate and the I7, that proves that the K = 4 state is really 
the 27 was given by means of 7 ray induced Doppler broadening following neutron capture 
by Borner et al. and using Coulomb excitation by Oshima et al. 12] and Hartlein et 

al. nj. 

fin 

The 27 states of very similar character were predicted |£ . |14j and observed in nearby nu- 
clides, 166 Er usin g Co ulomb excitation by Fahlander et al. 15| and using a (n, n/7) reaction 
by Garrett et al. 16], and 164 Dy using thermal neutron capture by Corminboeuf et al. ~v\ . 
In particular, the K = 27 state was also observed in 166 Er Sun et al. [3] studied the 
27 bands in 166 > 168 Er in terms of the triaxial projected shell model. This model gives the 
K = 4 2d states in a kinematical manner similar to Davydov-Filippov's asymmetric rotor 
n other mass regions, harmonic 27 states were observed in 232 Th by Korten 
ing Coulomb excitation and 106 ' 104 Mo by Guessous et al. [22|, |23j using a 
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et al. 
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spontaneous fission. The latter was discussed from a viewpoint of the X{5) symmetry 24]. 
These very limited number of observations indicate that it depends strongly on the under- 
lying single particle level structure whether the 27 states exist or not. This fact suggests 
that microscopic description of the collective vibrational excitations is mandatory. 

In odd- A nuclei, some numerical prediction for rare earth nuclides were made by Durand 



and Piepenbring 25j in terms of the multiphonon method prior to experimental observation 
The first experimental observation was made ten years later in 105 Mo by Ding et al. 2(| 
using a spontaneous fission. In these fission fragments, 104-106 Mo, rotational band members 
built on the 27 states were also populated. Soon after this, similar rotational bands were 
observed also in 103 Nb by Wang et al. 27] and 107 Tc by Long et al. j^ . These nuclei 
exhibit anharmonicity -S27/-E17 < 2 in 105 Mo and 103 Nb, while > 2 in 107 Tc, where level 
energies are measured from the corresponding zero phonon states. These are the 27 states 
observed so far. The first theoretical calculation to the observed 27 band in the odd- A 
nucleus, 103 Nb, was done by Sheikh et al. 29] in terms of the triaxial projected shell model. 
It was discussed there that the observed anharmonicity is difficult to be reproduced with 
the triaxial parameter that gives a good description of the I7 band. 

In the present paper, we take a complementary approach to the 27 band in 103 Nb; the 
particle vibration coupling (PVC) calculation based on the RPA phonons constructed in the 
rotating frame allowing possible static triaxial deformation. 
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II. THE MODEL 



In the model adopted in this study, elementary modes are quasiparticles and 7 vibrational 
RPA phonons excited on top of common vacuum as in traditional calculations, however, the 
vacuum is the yrast configuration of an even-even nucleus rotating with a frequency u TOt . 
This vacuum configuration can be either zero quasiparticle, two quasiparticle and so on, 
seen from the non-rotating ground state. Excitations are labeled by the signature quantum 
number r = exp (— iira), I = a + even, appropriate for rotating reflection-symmetric objects. 
Two kinds of 7 vibrational excitations exist, denoted by xK ± s with r = ±1, respectively. 
The phonon space is limited to zero, one and two phonon states. The vacuum mean field 
is rotating and static triaxial deformation is also possible; this means that the f2-mixing in 
quasiparticles and the K- mixing in RPA phonons are naturally taken into account. Here Q 
is the projection of the single particle angular momentum to the third axis. 

The formulation is summarized as follows. We begin with a one-body Hamiltonian in the 
rotating frame, 

ti = h- u mt J x , (1) 

h = h m -A T {P} + P T )-X T N T , (2) 

hm = ^ + \M(u%? + + uy) + v ls \ ■ s + v n {\ 2 - (\ 2 } No J. (3) 

In Eq. (j2J), P T is the pair annihilation operator, r = 1 and 2 denote neutron and proton, 
respectively, and the chemical potentials A r are determined so as to give the correct average 
particle numbers (N T ). The oscillator frequencies in Eq. ([3]) are related to the quadrupole 
deformation parameters €2 and 7 in the usual way. They, along with the pairing gaps 
A T , are determined from experimental information. The orbital angular momentum 1 in 
Eq. (jHJ) is defined in the singly stretched coordinates x' k = y~^ x k and the corresponding 
momenta, with k = 1, 2 and 3 denoting x, y and z, respectively. We apply the RPA to the 
residual pairing plus doubly stretched quadrupole-quadrupole (Q" ■ Q") interaction between 
quasiparticles. It is given by 

H mt = - £ g t pip t - \ £ ^qT'qT - \ E - { k ] qT )] qT\ (4) 

t=1,2 X=0,l,2 X=l,2 

where the doubly stretched quadrupole operators are defined by 

Q"k = Qx{xk -> x'u = —Xk), (5) 



and those with good signature are 



^2(1 + ^x0) 



(Q K ± Q^ K ) , 



(6) 



and P T is defined by subtracting the vacuum expectation value from P T . Among RPA modes 
determined by the equation of motion, 



[ti + H mt , Xl] RpA = u n Xl, 



(7) 



we choose the 7 vibrational phonons, n = 7(±), which have outstandingly large K = 2 
transition amplitudes 



\T, 



(±)l 



K 



n (±) x ] 



K »^ 7 (±) 



(8) 



(9) 



The particle vibration coupling Hamiltonian takes the form 

#couplc(7) = E^W^ ( X 7(+) a i°^ + X 7(+)4°m) 
+ J^ A 7(+)(^) (^(+)4^ + ^7(+)4%) 

+ E A 7(-)(^) (^ 7 (-)4<V + ^ 7 (-) a J t a ^) ' 

where and /i denote quasiparticles with r — —i and +i, respectively. The coupling vertices 
are given by 

a 7(+) (H = - E 4 +) ^ (+) q2 + V), 

K=0,l,2 

(10) 

A 7( „)M = - E 

K=l,2 

x"=i,2 

where Q"^ denotes quasiparticle scattering matrix elements that do not contribute to 



5 



RPA phonons. Eigenstates of the Hamiltonian thus specified at each u Iot take the form 

ixi> = E^V) 4W 
it 

+ X>j 3 Vr) <xjw + E>j 3) (£?) 4^ 



r\'\,^)—alx\X^\ 



(11) 



for the r 



-i sector, 



where 7 and 7 abbreviate 7(+) and 7(— ), respectively, and \<f>) is the rotating vacuum 
configuration. Those for the r = +i sector take a form similar to above, except that the 
suffices /i are to be replaced by p,. 



This model was first developed for studying the signature dependence o: 



and E2, Ml transition rates in one quasiparticle (zero phonon) bands 30 



the level energies 



applied to study the E2 intensity relation [4j, |32] in the I7 rotational bands 33j. The present 



31], and then 



study is the first application to the 27 states in rotating odd- A nuclei. By construction of 
the model space, this model is applicable up to the energy region < E7 lqp + 2A; otherwise 
non-collective 3qp states dominate over the collective states. 



III. RESULTS AND DISCUSSIONS 



Three new rotational bands that feed the 7r[422 5/2 + ] ground band of 4i 3 Nbgo were ob- 
served by a very recent 7 ray study of spontaneous fission fragments from 252 Cf |27j. Among 
them, the one that is built on the 9/2 + bandhead was assigned to the K = Q + 2 sequence 
of the I7 bands. One of the others that is built on the 13/2 + bandhead was assigned to 
the K = Q + 4 sequence of the 27 bands. This is the second observation of the 27 band in 



odd-A nuclei, and to which the first theoretical calculation [29] was reported. In the present 
cranked shell model calculation, diagonalization is performed in the five major shells, A osc = 
2-6 for the neutron and 1-5 for the proton with the Nilsson parameters vi s and vu taken 
from Ref. 3_4j. The three mean field parameters, which are assumed to be u rot independent 
for simplicity, the pairing gaps A n = 1.05 MeV, = 0.85 MeV and the deformation e-i 



0.31 are adopted from the experimental analyses [22|, |27j. The triaxiality 7 is chosen so as to 



reproduce the measured signature splitting of the ground band in the PVC calculation. The 
chosen value, 7 = —7°, gives overall reproduction of the signature splitting in the rotating 
frame aside from near the bandhead as shown in Fig. Q](a). Note here that it was reported 
that the signature splitting can be reproduced without invoking 7 deformation in the an- 
cestral model 35| of the triaxial projected shell model. This sug gest s that the appropriate 



value of the 7 deformation is model dependent as noticed in Ref. [36| and discussed below. 

The strengths of the residual doubly stretched quadrupole interaction are determined as 
follows; in the reference configuration with u Iot = and 7 = in which K is a good quantum 
number, determined to reproduce within the RPA the observed 7 vibrational 

energy in the adjacent 104 Mo, 0.812 MeV. If fully collective (3 vibration exists, can be 
determined to reproduce its energy but since the collective character of the observed 0^ 
is not clear, is set equal to And k[ +) = k[- ] is determined so as to make the 

energy of the Nambu-Goldstone mode zero. Those of the residual pairing interaction are 
determined to reproduce the adopted pairing gaps. Then the RPA calculation is performed 
with 7 = —7°. The obtained u rot dependence and signature splitting 0, Q of the excitation 
energy of 7 vibration is weak as shown in Fig. [T^b). 

Using these quantities, -ff C ou P ie(7) is diagonalized in the space of dimension 15 (number 



of quasiparticle states with N osc = 4) x 6 (lqp, lqp<8>7(+), lqp <g> 7(— ), lqpC8>7(+) 



7(+), lqpCg>7(— ) £g> 7(— ), lqp <8> 7(+) <S> 7(— ), with bar denoting the opposite signature) 
for each signature sector. Distribution of the strength (probability in the wave function) of 
the 7r[422 5/2 + ] <g> 7 vibration(s) that shows the collectiveness of each eigenstate directly is 
presented in Fig. [2] for the case of the favored signature r — —i. The result for the unfavored 
r = +i is similar. 

In the favored (r = —i) sector, two dominantly I7 eigenstates are obtained as a result of 
the interaction between the /®7(+) and the u® r y{— ) basis states, hereafter / and u denote 
the favored and unfavored lqp states originating from the 7r[422 5/2 + ], respectively. There 
is no general rule of the correspondence between these two bands in the signature scheme 
and the K = Q ± 2 bands in the K scheme. However, since states with the lower K have 
lower intrinsic energies than those with higher K and the same /, the obtained lower band 
can be identified with the K = Q — 2 band. Actually, in the previous analysis of the E2 



intensity relation 



33|, correspondence between the observed and the calculated states was 



established in this manner. The present result indicates that collectivity does not fragment 
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FIG. 1: (Color online) (a) Experimental and calculated signature splitting in the 7r[422 5/2 + ] one 
quasiparticle band as functions of the rotation frequency. Theoretical curve is the result of the 
particle vibration coupling rather than the cranked shell model alone. The latter (not shown) is 
larger than the former, (b) Excitation energies of the two types of 7 vibrational RPA phonons as 
functions of the rotation frequency. 



much and two bands are almost parallel as in the 165 Ho case 33[. But the higher band is 
slightly more collective and purer. This is consistent with the observation at u vot = that 
states with lower K are affected more by interaction with other states 25|, |38J. Actually 
only the K = Q + 2 sequence was observed in many cases. 
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FIG. 2: (Color online) Distribution of the strength (probability in the wave function) of the lqp, 

lqp(g>7 and lqp(8>7 (g) 7 components in the favored signature (r = — i), \ip^(fi)\ 2 (blue dotted), 

\i/A 3 \[ij)\ 2 + IV^HAt)! 2 (green dashed) and |V> (5) (/n7)| 2 + IV> (5) (M77)| 2 + |V> {5) (^77)1 2 (red solid), 

respectively, at various rotation frequencies, w ro t = - 0.35 MeV. Here lqp (fi) means the one 

quasiparticle states originating from the 7r[422 5/2 + ] at oj rot = 0. 

9 



Fragmentation of the 27 components is expected to depend sensitively on how other 
quasiparticle states distribute. In the present calculation, the 27 strength concentrates 
mainly on three eigenstates that locate at < -Ei qp + 2A P . They are obtained as a result of 
the interaction among the / (g) 7(+) <8> 7(+), the / ® 7(— ) <S> 7(— ) and the u (8> 7(+) ® 7(— ) 
basis states. With the same rule as for the I7 case, they can be identified with K — — 
Q and Q + 4 from the lower. At u vot = the higher K states are the more collective as 
expected. But small rotation immediately mixes their collectivity. Then three collective 
bands run keeping almost the same collectivity. 

In order to see how these collective states interact with other non-collective states, all the 
eigenstates of H conp \ e {^) that locate lower than 3 MeV are shown in Fig. EJ^a). At u TOt = 0, 
there are 2nd - 4th iV osc = 4 lqp states between the 17s and 27s. These are expected to 
correspond to some of the levels that were known to be populated by the 103 Zr f3 decay but 
spin and parity have not been assigned 39J. One of them interacts with the 27s depending 
on w ro t. Aside from this, the 27s keep their collective characters as shown in Fig. |2j 

Next the obtained collective states are compared with observed ones in Fig. ]3](b). This 
figure shows that the lqp and the I7 (K = Q + 2) are reproduced almost perfectly. On the 
other hand, the calculated 27s obviously locate higher than the observed one. The most 
probable reason of this is ignoring higher lying states, that is, their effect to push down 
the 27 states. Although one may be afraid that inclusion of more states would lead to 
fragmentation of collectivity, we expect this not to occur since the number of nearby states 
is still small. 

Here we discuss the adopted triaxial deformation 7 in Ref. 29[ and the present calculation. 
It has two physical aspects, its sign and its absolute value. The former is originally the 
notion in the rotating mean field model. Although the projection model in its original 
framework does not distinguish the sign of 7 to our knowledge, Ref. |36| proposed a method 
to find the main rotation axis. This made it possible to relate the result of the projection 
calculation to the sign of 7. In Ref. [29 1. the adopted e = 0.3 and e' = 0.16 correspond 
to 7 = ±28°. Although its sign is not described, odd- A nuclei with an Q = 5/2 odd 
nucleon show 7 < rotation due to its shape driving effect in addition to the approximately 
irrotational property of the even-even core in general. Then we regard the sign is consistent 
with ours. As for the latter, its absolute value, the origin of the difference is that of the 
generating mechanism of the 7 band itself as pointed out in Ref. [181 ] . It is constructed as 
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a vibrational phonon excitation built on top of the yrast configuration in our model. In 
contrast, in the projection model, it is generated as a rotational excitation in the sense of 
the asymmetric rotor. Consequently appropriate 7 deformation is not necessarily the same. 
The latter calls for a larger |7|; for example, I7I = 25° was adopted for 168 Er [18], which 
is regarded as a typical axially symmetric nucleus |4J. Note here that large fluctuations in 
the 7 direction were shown for nuclei of that class, for example, in Ref. j40j. These results 
indicate that the model dependence of the adopted value of 7 deformation is not unphysical. 
It is an interesting subject to look into more this longstanding problem of how these two, 
vibrational/rotational, descriptions are related, but it is beyond the scope of the present 
paper. 

Last but not least, we mention another type of collective motion for which two phonon 
excitation was observed, the wobbling motion predicted by Bohr and Mottelson [ 4| , Janssen 



and Mikhailov 411 and Marshalek 



42j . and first observed by 0degard et al. 



43|. This 



is a small amplitude fluctuation of the rotation axis of triaxially deformed nuclei and has 
the same quantum number, r = —1, as the odd spin members of the 7 vibrational bands. 
Two phonon wobbling bands were observed in 163 Lu by Jensen et al. 44j and 165 Lu by 



Schonwafier et al. 



45]. Their excitation energies in the rotating frame exhibit anharmonicity 



such that the two phonon states locate lower than the twice the energies of the one phonon 



states. The present author and Ohtsubo 



46] argued that this anharmonicity indicates the 



softening of the collective potential surface, that is, the precursor of the second order phase 
transition from the principal axis rotation to the tilted axis rotation. This is completely 
parallel with the situation that the softening of the 7 vibration is the precursor of the second 
order phase transition from the axial symmetric deformation to static triaxial deformation. 
The new vacuum after the phase transition can accommodate anharmonic vibration either 
i?2- P honon/-Ei-phonon > 2 or < 2. Compare the potential energy surfaces: Fig.3 in the first 
item of Ref. [8] and Fig.3 in Ref. 46 1; in general, stiff potentials as the former lead to high 
-£'2-phonon/-E'i -phonon ratios while soft ones as the latter do low ratios. 

To summarize, two phonon 7 vibrational bands in rotating odd- A nuclei were observed 



recently in three nuclides 



26j-|28| . After these observation, the first theoretical calculation by 



. In the present 



means of the triaxial projected shell model to one of them was reported 
paper, we have applied to the same nuclide a completely different, mean field based model, 
which was previously developed for the signature dependent properties of one quasiparticle 
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(zero phonon) bands and later utilized for the E2 intensity relations in one phonon bands. 
In the present calculation, after confirming almost perfect reproduction of the zero and one 
phonon states, we concentrated on the effect of the odd quasiparticle on the distribution of 
the two phonon collectivity. The obtained K = Q + 4 state is almost pure at u rot = 0, but 
small rotation immediately delivers the strength to two other sequences; consequently three 
collective sequences keep about 60% collectivity up to higher spins. This indicates future 
experimental observation of more than one sequence of the 27 bands as well as the other 
sequence of the I7. The restriction of the model space mentioned above would be the reason 
of the result that the calculated 27 states locate higher than observed. 

There are some to improve the present calculation: First of all, the model space should 
be enlarged to more than two phonon states. Further, extension of the microscopic approach 
by means of the selfconsistent collective coordinate method to rotating even and odd- A 
systems is promising in future. 
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FIG. 3: (Color online) (a) Calculated eigenstates of -ff C oupie(7) at each rotation frequency. Cor- 
responding to Fig. [2j those with large (> 40%) 7r[422 5/2 + ] <g> 0, 1,27 components are marked by 
symbols. At u; rot = 0.25, 0.275 and 0.325 MeV, accidental fragmentation caused by interaction 
with a positive sloping lqp state occurs, (b) Similar to (a) but the dominantly 7r[422 5/2 + ] ® 0, 1, 27 
states are compared with experimental data (curves) converted to the rotating frame by using the 



Harris parameters J = 15.45 MeV" 1 and Ji = 81.23 MeV" 3 that fit the yrast band of 104 Mo 



23(]. 
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